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Abstract 

The circular Radon transform integrates a function over the set of 
all spheres with a given set of centers. The problem of injectivity of 
this transform (as well as inversion formulas, range descriptions, etc.) 
arises in many fields from approximation theory to integral geometry, 
to inverse problems for PDEs, and recently to newly developing types 
of tomography. A major breakthrough in the 2D case was made sev- 
eral years ago in a work by M. Agranovsky and E. T. Quinto. Their 
techniques involved microlocal analysis and known geometric proper- 
ties of zeros of harmonic polynomials in the plane. Since then there 
has been an active search for alternative methods, especially the ones 
based on simple PDE techniques, which would be less restrictive in 
more general situations. The article provides some new results that 
one can obtain by methods that essentially involve only the finite speed 
of propagation and domain dependence for the wave equation. 



1 Introduction 

The circular Radon transform integrates a function over the set of all spheres 
with a given set of centers. Such transforms have been studied over the years 
in relation to many problems of approximation theory, integral geometry, 
PDEs, sonar and radar imaging, and other applications |9] [T8]. 
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[ini-[ni, En|-Eni, EZI-ESI, EZI-lini)- Although significant progress has 
been achieved, some related analytic problems have proven to be rather hard 
and remain unresolved till now. A new wave of interest to such transform 
has been sparked by the recent development of the Thermoacoustic Tomog- 
raphy (TAT or TCT) as one of the promising methods of medical imaging 
(e.g., |2I|5lSZ|-|in])- The TAT procedure can be described as follows: a short 
microwave or radiofrequency electromagnetic pulse is sent through the bio- 
logical object. At each internal location x certain energy H{x) is absorbed. 
The cancerous cells happen to absorb several times more MW (or RF) energy 
than the normal ones, which means that significant increases of the values 
of H{x) are expected at tumorous locations. The absorbed energy, due to 
resulting heating, causes thermoelastic expansion, which in turn creates a 
pressure wave. This wave can be detected by ultrasound transducers placed 
outside the object. It has been shown that here one effectively measures 
the integrals of H{x) over all spheres centered at the transducers' locations. 
In other words, one needs to invert the mentioned above generalized Radon 
transform of H ("generalized," since integration is done over spheres). It is 
clear from the dimension considerations that it should be sufficient to run 
the transducers along a curve in the case of a 2D problem or a surface in 
3D. The most popular geometries of these surfaces (curves) that have been 
implemented are spheres, planes, and cylinders 

The central problems that arise in these studies are: uniqueness of re- 
construction, reconstruction formulas and algorithms, stability of the recon- 
struction, description of the range of the transform and incomplete data 
problems. 

All these questions have been essentially answered for the classical Radon 
transform that arises in X-ray CT, Positron Emission Tomography (PET), 
and Magnetic Resonance Imaging (MRI) fl^ 177] . However, they are much 
more complex and not that well understood for the circular Radon transform 
that arises in TAT. 

This paper contains some new approaches and results concerning the 
uniqueness problem. The reader should be aware that for the currently prac- 
tically used geometries of TAT the uniqueness issue has been resolved. E.g., 
for the spherical location of the centers (transducers) uniqueness follows for 
instance from Corollary first proven in [22^ (see also [21 E] and references 
therein). For the planar location, it has been known for a long time [201 IH] 
that only odd functions with respect to this plane cannot be reconstructed 
from the spherical integral data. However, the complete understanding of the 
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uniqueness problem for general locations of the transducers remains elusive 
(especially in dimensions higher than two). The aim of this paper is to make 
progress in filling this gap by obtaining new uniqueness results, as well as by 
reproving some known results by simpler means, which makes them easier to 
extend to higher dimensions and other geometries. 

The results of this paper were presented at the special sessions on to- 
mography at the AMS Meetings in Binghamton, NY in October 2003 and 
in Lawrenceville, NJ in April 2004 and at the Inverse problems workshop at 
IPAM in November 2003. 

The next section contains the mathematical formulation of the problem 
and its brief history. The following section contains the main results of this 
paper. It is followed by sections containing further remarks and acknowl- 
edgements. 

2 Formulation of the problem and known re- 
sults 

The discussion of the previous section motivates the study of the following 
"circular" Radon transform. Let f{x) be a continuous function on R", n > 2. 

Definition 1 The circular Radon transform of f is defined as 

Rf{p,r)= [ f{y)da{y), 

J\y~p\=r 

where da{y) is the surface area on the sphere \y — p\ = r centered at p E MJ^. 

In this definition we do not restrict the set of centers p or radii r. It 
is clear, however, that this mapping is overdetermined, since the dimension 
of pairs (p, r) is n + 1, while the function / depends on n variables only. 
This suggests to restrict the set of centers to a set (hypersurface) 5" C R", 
while not imposing any restrictions on the radii. We denote this restricted 
transform by Rs'- 

Rsf{p,r) = Rf{p,r)\,^s. 

Definition 2 The transform R is said to be injective on a set S (S is a 
set of injectivity^ if for any f G Cc(R"') the condition Rf(p,r) — for all 
r e R and all p E S implies / = 0. 
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In other words, S is a set of injectivity, if the mapping Rs is infective on 

Ce(M"). 

Here we use the standard notation Cc(M") for the space of compactly sup- 
ported continuous functions on M". The situation can be significantly differ- 
ent and harder to study without compactness of support (or at least some 
decay) condition [21 E|- Fortunately, tomographic problems usually yield 
compactly supported functions. 
One now arrives to the 

Problem 3 Describe all sets of injectivity for the circular Radon transform 
R on Cc(M"). 

This problem has been around in different guises for quite a while j3 El 
1231 1211 ■ The paper [3] contains a survey of some other problems that lead to 
the injectivity question for Rs- 

The first important observations concerning non-injectivity sets were made 
by V. Lin and A. Pincus [231 12^ and by N. Zobin [IT] . Their results imply 
in particular that if R is not injective on S", then S is contained in the zero 
set of a harmonic polynomial. Therefore we get a sufficient condition for 
injectivity: 

Corollary 4 Any set S <Z MJ^ of uniqueness for the harmonic polynomials 
is a set of injectivity for the transform R. 

In particular, this implies 

Corollary 5 If U dW^ is a hounded domain, then S = dU is a injectivity 
set of R. 

We will see later a different proof of this fact that does not use harmonicity. 

So, what are possible non-injectivity sets? Any hyperplane S is such a set. 
Indeed, for any function / that is odd with respect to S, one gets Rsf = 0. 
There are other options as well. In order to describe them in 2D, let us first 
introduce the following definition. 

Definition 6 For any N denote hyT.^^ the Coxeter system of N lines 
Lq, . . . , Ln-i in the plane^ : 

Lk = - cx)< t < oo}. 

^In the formula below we identify the plane with the complex plane C. 
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In other words, S^v is a "cross" of lines passing through the origin and 
forming equal angles t[/N . It is rather easy to construct a non-zero compactly 
supported function that is simultaneously odd with respect to all lines of a 
given Coxeter set. Hence, Sat is a non-injectivity set as well. Applying 
any rigid motion cu, one preserves non-injectivity property. It has been also 
discovered that one can add any finite set F preserving non-injectivity. Thus, 
all sets wEtv U F are non-injectivity sets. It was conjectured by V. Lin 
and A. Pincus that these are the only non-injectivity sets for compactly 
supported functions on the plane. Proving this conjecture, M. Agranovsky 
and E. Quinto established the following result: 

Theorem 7 '^31 The following condition is necessary and sufficient for a set 
S (ZM? to he a set of injectivity for the circular Radon transform on Cc(M^).- 
S is not contained in any set of the form ui{T,n) [JF, where u is a rigid 
motion in the plane and F is a finite set. 

The (unproven) conjecture below describes non-injectivity sets in higher 
dimensions. 

Conjecture 8 ^ The following condition is necessary and sufficient for S 
to be a set of injectivity for the circular Radon transform on Cc(]R").- 

S is not contained in any set of the form a;(E) [jF, where u is a rigid 
motion of M", S is the zero set of a homogeneous harmonic polynomial, 
and F is an algebraic subset in M" of co-dimension at least 2. 

The reader notices that for n = 2 this boils down to Theorem [7| 

The beautiful proof of Theorem [7| by M. Agranovsky and E. Quinto is 
built upon the following tools: microlocal analysis (Fourier Integral Opera- 
tors technique) that guarantees existence of certain analytic wave front sets 
at the boundary of the support of a function located on one side of a smooth 
surface (Theorem 8.5.6 in [19^), and known geometric structure of level sets 
of harmonic polynomials in 2D (e.g., [14j). These methods, unfortunately, 
restrict wider applicability of the proof. The microlocal tool works at an 
edge of the support and hence is not applicable for non-compactly-supported 
functions. On the other hand, the geometry of level sets of harmonic poly- 
nomials does not work well in dimensions higher than 2 or on more general 
Riemannian manifolds (e.g., on the hyperbolic plane). Thus, the quest has 
been active for alternative approaches since [3] has appeared. 
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It is instructive to look at alternative reformulations of the problem 
(which there are plenty |3j). There is a revealing reformulation [211221 that 
stems from known relations between spherical integrals and the wave equa- 
tion (e.g., [3 (201)- Namely, consider the initial value problem for the wave 
equation in M": 

utt - Au = 0, xeW,teR (1) 



u\t=o = 0, Ut\t=o = /• (2) 

Then 

^a;,t)= ^^_^^,^^^_^ y^ r(t2-0("-3)/2(i?/)(x,r)rfr, t > 0. 

Hence, it is not hard to show [3 that the original problem is equivalent to 
the problem of recovering ut{x,0) from the value of u{x,t) on subsets of 
5* X (— oo, oo). 

Lemma 9 A set S is a non-injectivity set for Cc(M") if and only if 
there exists a non-zero compactly supported continuous function f such that 
the solution u{x, t) of the problem ^\)-^ vanishes for any x E S and any t. 

Hence, non-injectivity sets are exactly the nodal sets of oscillating free infinite 
membranes. In other words, injectivity sets are those that observing the 
motion of the membrane over S gives complete information about the motion 
of the whole membrane. 

One can now try to understand the geometry of non-injectivity sets in 
terms of wave propagation. The first example of such a consideration was the 
original proof 122jj of Corollary El that did not use harmonicity (not known at 
the time). Let S = dU be a non-injectivity (and hence nodal for wave equa- 
tion) set, where ?7 is a bounded domain. Then on one hand, the membrane 
is free and hence the energy of the initial compactly supported perturbation 
must move away. Thus, its portion inside U should decay to zero. On the 
other hand, one can think that is a fixed boundary and hence the energy 
inside must stay constant. This contradiction allows one to conclude that 
in fact / = 0. The same PDE idea, with many more technical details, was 
implemented in [2: to prove the following statement: 

Theorem 10 l^IfUisa bounded domain in M", then S = dU is an injec- 
tivity set for R in the space L^(M'^) if q < 2n/ {n — 1). This statement fails 
when q > 2n/{n — 1), in which case spheres fail to be injectivity sets. 
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In spite of these limited results, it still had remained unclear what dis- 
tinguishes in terms of wave propagation the "bad" flat lines S in Theorem [7| 
that can be nodal for all times, from any truly curved S that according to 
this theorem cannot stay nodal. An approach to this question was found in 
the recent paper ^3] by D. Finch, Rakesh, and S. Patch, where in particular 
some parts of the injectivity results due to [21 were re-proven by simple PDE 
means without using microlocal tools and harmonicity: 

Theorem 11 llcif Let D be a bounded, open, subset o/M", n > 2, with a 
strictly convex smooth boundary S . Let T be any relatively open subset of 
S. If f is a smooth function on supported in D, u is the solution of the 
initial value problem (1), (2) and u{p, t) = for all t and p G F, then / = 0. 

Although this theorem follows from microlocal results in jS]^, its signifi- 
cance lies in the proof provided in JHl (that paper contains important results 
concerning inversion as well, which we do not touch here). 

The following two standard statements concerning the unique continua- 
tion and finite speed of propagation for the wave equation were the basis of 
the proof of the Theorem ITD in They will be relevant for our purpose 
as well. 

Proposition 12 11 ^Jf Let BJp) = {x G | \x—p\ < e}. Ifu is a distribution 
and satisfies ^ and u is zero on B^{p) x (— T, T) for some e > 0, andp G M", 
then u is zero on 

{(x, t) : |x — p| + \t\ < T}, 

and in particular on 

{(x, 0) : |x — p| < T}. 

Let now D be a bounded, open subset of M" with the boundary S. For 
points p, q outside let d{p^ q) denote the infimum of the lengths of all the 
piecewise paths in R" \ D joining p to q. Then d{p, q) is a metric on 
MP\D. For any point p 'm.W^\D and any positive number r, define Er{p) to 
be the ball of radius r and center at p in R" \ D with respect to this metric, 
i.e. 

Er{p) = {x eW\D : d{x,p) < r}. 

^Results of j3j make the situation described in Theorem ll ll impossible. since the support 
of / lies on one side of a tangent plane to T. See also Theorem [T^ and |25|. 
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Proposition 13 11 ^Jf Suppose D is a bounded, open, connected subset o/R", 
with a smooth boundary S . Let u be a smooth solution of the exterior problem 

utt- Au = 0, X eR"\D, teR 
u = h on S X R. 

Suppose p is not in D, and to < ti are real numbers. If u{.,to) and ut{.,to) 
are zero on Et^^toip) o-nd h is zero on 

{(x, t) : X & S,tQ < t < ti, d{x,p) < ti — t}, 

then u{p, t) and Ut{p, t) are zero for all t E [to, ^i)- 

3 Further injectivity results by PDE means 

We will show now how simple tools similar to the Propositions ^1 and 
Uni namely finite speed of propagation and domain of dependence for the 
wave equation allow one to obtain more results concerning geometry of non- 
injectivity sets, as well as to re-prove some known results with much simpler 
means. The final goals were to recover the full result of j^] in 2D and to prove 
its analogs in higher dimensions and for other geometries (e.g., hyperbolic 
one) using these simple means. Albeit this goal has not been completely 
achieved yet, we can report some progress in all these directions. 

Let us start with some initial remarks that will narrow the cases we need 
to consider. First of all, one can assume functions / as smooth as we wish, 
since convolution with smooth radial moUifiers does not change the fact that 
Rsf = (e.g., Secondly, according to the results mentioned before, any 
non- injectivity set S in the class of compactly supported functions is con- 
tained in an algebraic surface that is also a non-injectivity set. It is rather 
straightforward to show that the same is true for functions that decay expo- 
nentially. Thus, considering only exponentially decaying functions, 
one does not restrict generality by assuming from the start alge- 
braicity of S. It is known [1^ that algebraic surfaces of co-dimension higher 
than 1 are automatically non-injectivity sets. Thus, we can restrict our atten- 
tion to algebraic hypersurfaces 5* of only. Any set that is not algebraic (or 
rather, is not a part of such an algebraic surface) is automatically an injec- 
tivity set. So, when trying to obtain necessary conditions for non-injectivity, 
confining ourselves to the case of algebraic hypersurfaces solely we do not 
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lose any generality. One can also assume irreducibility of that surface, if this 
helps. When needed, one can also exclude the case of closed hypersurfaces, 
since according to Corollary El those are all injectivity sets. 

Our goal now is to exclude some pairs (5,/), where S is an algebraic 
surface and / is a non-zero function as possible candidates for satisfying the 
non- injectivity condition Rsf = 0. We will do this in terms of geometry 
of the support of function /. Notice that Theorem ^2 does exactly that 
when S contains an open part of the boundary of a smooth strictly convex 
domain where / is supported. Theorem [7J on the other hand excludes all 
compactly supported /'s in R^, unless S = cuSat. Similarly, Theorem ITUl 
excludes boundaries S of bounded domains when / is in an appropriate 
space Lp{W). 

Let S be an algebraic hypersurface (which can be assumed to be irre- 
ducible if needed) that splits into connected parts , j = 1, ...,m. One 
can define the interior metric in as follows: 

d^(p,q) = m/{length of 7}, (3) 

where the infimum is taken over all C^-curves 7 in joining points p,q ^ HK 

Theorem 14 Let S and be as above and f G C(]R") be such that Rsf = 
0. Let also x G , where is the closure of . Then 

dist{x, supp f n H^) = disP{x, supp f fl H^) , . 

< dist{x, supp f n H'^), k j, ^ ' 

where distances disP are computed with respect to the metrics d^ , while dist 
is computed with respect to the Euclidean metric in MJ^. 
In particular, for x G S* and any j 

dist{x, supp f n H^) = disP{x, supp f fl H^) = dist{x, supp /). (5) 

Thus, the expressions in ^ in fact do not depend on j = 1, ...,m. 

Remark 15 Notice that under the condition of algebraicity of S the theorem 
does not require the function f to be compactly supported and in fact imposes 
no condition on behavior of f at infinity. On the other hand, as it has been 
mentioned before, if f decays exponentially, then the algebraicity assumption 
does not restrict the generality of consideration. 
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Proof of the theorem. Notice first of all, that the function d^{p,x) has 
gradient \Vxd^ {p,x)\ < 1 a.e.^ 
Let us prove now the equality 

dist{x, supp f n H^) = dist^{x, supp f r\ W). (6) 

Since d^{p,q) > \p — q\, it is sufficient to prove that the left hand side ex- 
pression cannot be strictly smaller than the one on the right. Assume the 
opposite, that 

dist{x, supp f n H^) = di < d2 = dist^ {x, supp f H H^). (7) 

Pick a smaller segment [d^, d^] C ((ii, ^2)- Then, by continuity, for any point 
p in a small ball B C near x (not necessarily containing x, for instance 
when X ^ S) one has 

dist{p, supp / n H^) <d'i<di< dist\p, supp / fl H^). (8) 

For such a point p, consider the volume V in the space-time region x M 
bounded by the space-like surfaces Si given by t = and S2 described as t = 
(f){x) = T — d^{p,x), and the "vertical" boundary SxR. Here r < {ds + d4)/2. 
Consider the solution u{x,t) of the wave equation problem (IH)-(|11) with the 
initial velocity /. Then, by construction, this solution and its time derivative 
are equal to zero at the lower boundary t = and on the lateral boundary 
5 X R. Hence, by the standard energy computation (integrating the equality 
uDu = 0, see, e.g.. Section 2.7, Ch. 1 in 0) we conclude that m = in V^. 
For the reader's convenience, let us provide brief details of the corresponding 
calculations from jS]: Since □« = 0, u = = on Ei, and u\s = for all 
times, we get by integration by parts 

= JutDudxdt= J ^{\Vu\^ + u'^ + 2ut'V(f)-Vu)dx 
= i / {\V{u{xMx))\' + {l-\V<P\')ut{x,<P{x)Y)dx. 

(t>{x)>0 

Since |V(/)| < 1, we conclude that 

j (|V(m(x,0(x))P) rfx = 

<l>{x)>0 

^In order to justify legality of the calculation presented below, one can either use 
geometric measure theory tools, as in J^], or just notice that due to algebraicity of S, the 
function (p, x) is piece-wise analytic. 



10 



and hence u is constant on S2. Taking into the account the zero conditions 
on S and Si, one concludes that m = on E2, and hence in V . 

In particular, u{p,t) = for all p E B and \t\ < {d^ + d4)/2. Notice that 
(^3 + d4)/2 > (is. Now applying Proposition [T21 to the wave equation in the 
whole space, we conclude that 

dist{p, supp /) > ds, (10) 

and hence 

dist{p, supp / n H^) > ^3, (11) 
which is a contradiction. This proves ©. It is now sufficient to prove 

dist{x, supp / n H^) < dist{x, supp / n H'') (12) 

for k ^ j. This in fact is an immediate consequence of (jl(J|) . Alternatively, we 
can repeat the same consideration as above in a simplified version. Namely, 
suppose that 

dist{x, supp / n > d2> di> dist{x, supp / fl H^) (13) 

for a point x G fl S, and hence for all points p in a small ball in . 
Consider the volume V in the space-time region x M bounded by the 
space-like surfaces t = and t = d2 — \x—p\ {p fixed in the small ball) and the 
boundary S* x M. Consider the solution u{x,t) of the wave equation problem 
with the initial velocity /. Then, by construction, this solution and 
its time derivative are equal to zero at the lower boundary t = and on the 
lateral boundary S" x M. Hence, by the same standard domain of dependence 
argument (see, e.g.. Section 2.7, Ch. 1 in jH]) we conclude that m = in V^. In 
particular, u{p,t) = for all p G -B and \t\ < ^2- Now applying Proposition 
1121 to the wave equation in the whole space, we conclude that 

dist{p, supp /) > d2, 

and hence 

dist{p, supp / n if ^) > (i2, (14) 

which is a contradiction. □ 

We will now show several corollaries that can be extracted from Theorem 

M 
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Corollary 16 Let f be continuous and S CW^ be an algebraic hypersurface 
such that Rsf = 0. Let L be any hyperplane such that L fl supp / 7^ and 
such that supp f lies on one side of L. Let x G LCisupp f and Tx be the open 
ray starting at x, perpendicular to L, and going into the direction opposite to 
the support of f . Then either C S (and hence, the whole line containing 
Tx belongs to S), or does not intersect S. 

Proof. Assuming otherwise, let p & (1 S and be the connected com- 
ponents of M'^yS* such that p belongs to their closures. Since x is the only 
closest point to p in the support of /, Theorem implies that for any j 
there exist paths t^ joining x and p through and such that the length of 
tends to \x — p\ when e — 0. This means that these paths converge to 
the linear segment [x,p]. Hence, this segment belongs to for any j, and 

thus to nH\ which is a part of S. We conclude that the segment [x,p], and 
j 

then, due to algebraicity of S, the whole its line belongs to S. This proves 
the statement of the corollary. □ 

One notices that a similar proof establishes the following 

Corollary 17 Let f be continuous and S C M" be an algebraic hypersurface 
such that Rsf = 0. Suppose p E S is such that p does not belong to supp f 
and there exists unique point x in supp f closest to p. Then S contains the 
whole line passing through the points x and p. 

Let S CW"-. For any points p, g G — S* we define the distance ds{p, q) 
as the infimum of lengths of paths in M" — S connecting these points. 
Clearly ds{p, q) > \p — q\- Using this metric, we can define the corresponding 
distances dists from points to sets. 

Theorem 18 Let a set S C M" and a non-zero function f G C(M") ex- 
ponentially decaying at infinity be such that Rsf = 0. Then for any point 
peW'-S 

dists{p, suppf) = dist{p, suppf). (15) 

The same conclusion holds for any continuous function, if one assumes that 
S is an algebraic hypersurface. 

Proof. Assume that p5|) does not hold, i.e. 

(iist5(p, supp/) > (iist(p, supp/). 
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As it has been mentioned before, under the conditions of the theorem, we can 
assume S" to be a part of an algebraic surface S for which R-^f = 0. Let S 
divide the space into parts . Then, in notations of the previous theorem, 
we have 

dist^p, supp / n H^) > dists{p, supp/) (16) 

and hence 

disP{p, supp / n H^) > dist{p, supp/). (17) 

This, however, contradicts Theorem ITU □ 

Let us formulate another example of a geometric constraint on pairs S, f 
such that Rsf = 0^ 

Theorem 19 Let S C be a relatively open piece of a -hypersurface 
and f G Cc(M") be such that Rsf = 0. // there is a point po E S such that 
the support of f lies strictly on one side of the tangent plane Tp^^S to S at 
Po, then f = 0.^ 

Proof of the theorem. Let us denote by i^p(supp /) the convex cone with 
the vertex p consisting of all the rays starting at p and passing through the 
convex hull of the support of /. Then Kp^{supp /), due to the condition of 
the theorem, lies on one side of Tp^S, touching it only at the point Pq. Let us 
pull the point po to the other side of the tangent plane along the normal to 
a nearby position p. Then it is easy to see that for p sufficiently close to po, 
all rays of the cone Kp{supp f) will intersect S. This means in particular, 
that for this point p we have dists{p, supp f) > dist{p, supp f). According 
to Theorem ^[ this implies that / = 0. □ 

Corollary 20 Let S C M" be an algebraic hypersurface and f G Cc(]R"). 
// Rsf = 0, then every tangent plane to S intersects the convex hull of the 
support of f . 

The above results present significant restrictions on the geometry of the 
non-injectivity sets S and of the supports of functions / in the kernel of Rs- 
One can draw more specific conclusions about these sets. The statement be- 
low was proven in ^ by using the geometry of zeros of harmonic polynomials, 
which we avoid. 

similar statement in the case of analytic surfaces S was announced in n25j for 
distributions /. The proof is claimed to be based upon microlocal analysis. 
^This implies, in particular, Theorem II II 
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Proposition 21 Let S G M."^ be an algebraic curve such that Rsf = for 
some non-zero compactly supported continuous function f . Then S has no 
compact components, and each its component has asymptotes at infinity. 

Proof. CorollaryElexcludes bounded components. So, we can think that S is 
an irreducible unbounded algebraic curve. Existence of its asymptotes can be 
shown as follows. Let us take a point p G S and send it to one of the infinite 
ends of S. According to Corollary EDI every tangent line TpS intersects the 
convex hull of the support of /, which is a fixed compact in R^. This makes 
this set of lines on the plane compact. Hence, we can choose a sequence 
of points pj such that the lines Tp^S converge to a line T in the natural 
topology of the space of lines (e.g., one can use normal coordinates of lines 
to introduce such topology). This line T is in fact the required asymptote. 
Indeed, let us choose the coordinate system where T is the x-axis. Then the 
slopes of the sequence Tp. S converge to zero. Due to algebraicity, for a tail 
of this sequence, the convergence is monotonic, and in particular holds for 
all p e 5 far in the tail of S. Let us for instance assume that these slopes 
are negative. Then the tail of S is the graph of a monotonically decreasing 
positive function. This means that 5* has a horizontal asymptote. This 
asymptote must be the x-axis T, otherwise the intercepts of Tp.S would 
not converge to zero, which would contradict the convergence of Tp.S to T. 
□ 

The next statement proves the Agranovsky-Quinto Theorem [7| in some 
particular cases. In order to formulate it, we need to introduce the following 
condition: 

Condition A. Let i^' be a compact subset of M". We will say that the 
boundary of K satisfies condition A^, if there exists a positive number tq 
such that for any r < Tq and any point x in the infinite connected component 
of M" \ K such that dist (x, K) = r there exists a unique point k on K such 
that \x — k\ = r. 

Examples of such sets are convex sets (where tq > can be chosen arbi- 
trarily) and sets with a boundary (where ro should be sufficiently small). 

Theorem 22 Let S CM.'^ and /(^ 0) G ^^(M^) be such that Rsf = 0. If the 
external boundary of the support of f (i.e., the boundary of the infinite com- 

^This condition essentially restricts the curvature of the boundary from below. 
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ponent of the complement of the support) is connected and satisfies Condition 
A, then S C cuSjv U F in notations of Theorem^ 

The conditions of the theorem are satisfied for instance when the support 
of f contains the boundary of its convex hull, or when the support's external 
boundary is connected and of the class C^. 

Proof. First of all, up to a finite set, we can assume that S is an algebraic 
curve. Since the external boundary of the support is assumed to be con- 
nected. Theorem implies that any irreducible component of 5* must meet 
any neighborhood of the support of /. If we take the neighborhood of radius 
r < tq, then each point on S in this neighborhood will have a unique closest 
point on supp /. Applying now Corollary [T71 we conclude that 5* consists of 
straight lines Lj intersecting the support. It is known that any straight line L 
is a non-injectivity set, but the only functions annihilated by Rl are the ones 
odd with respect to L (e.g., [21 El 1201)- Hence, / is odd with respect to all 
lines Lj. In particular, every of these lines passes through the center of mass 
of the support of /. Hence, lines Lj form a "cross" It remains now to show 
that the angles between the lines are commensurate with vr. This can also be 
shown in several different ways. For instance, this follows immediately from 
existence of a harmonic polynomial vanishing on S. Another simple option 
is to notice that if this is not the case, then there is no non-zero function 
that is odd simultaneously with respect to all the lines. □ 

Exactly the same consideration as above shows that in higher dimensions 
the following statement holds: 

Proposition 23 Let 5 C M" and /(^ 0) G Cc(K") be such that Rsf = 0. 
// the external boundary of the support of f (i.e., the boundary of the infi- 
nite component of the complement of the support) is connected and satisfies 
Condition A, then S is ruled (a scroll)^. 

The conditions of the theorem are satisfied for instance when the support 
of f contains the boundary of its convex hull, or when the support's external 
boundary is connected and of the class C^. 

^One can prove that all these lines pass through a joint point also in a different manner. 
Indeed, due to oddness of /, each line is a symmetry axis for the support of /. Then, 
considering the group generated by reflections through these lines, one can easily conclude 
that if they did not pass through a joint point, then the support of / must have been 
non-compact. 

ruled surface, or a scroll is the union of a family of lines (e.g., {361 ) 
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Remark 24 // we could also show that all these lines pass through the same 
point, then this would immediately imply, as in the previous proof, the validity 
of Conjectured^ for this particular case. 

4 Additional remarks 

1. M. Agranovsky and E. T. Quinto have written besides [3], several other 
papers devoted to the problem considered here. They consider some 
partial cases (e.g., distributions / supported on a finite set) and varia- 
tions of the problem (e.g., in bounded domains rather than the whole 
space). See jTHUEHn] for details. 

2. One of our goals was to obtain the complete Theorem [71 the main 
result of [3] by simple PDE tools, avoiding using the geometry of zeros 
of harmonic polynomials and microfocal analysis (or at least one of 
those), as well as to prove its analogs in higher dimensions and for other 
geometries (e.g., hyperbolic one). Although we have not completely 
succeeded in this yet, the results presented (e.g.. Propositions 1711 and 
1221 and Theorem are moving in this direction. 

3. The PDE methods presented here in principle bear a potential for con- 
sidering non-compactly-supported functions. In order to achieve this, 
one needs to have qualitative versions of statements like Proposition 
Eland Theorem ITHl where instead of just noticing whether a wave has 
come to certain point at a certain moment (which was our only tool), 
one controls the amount of energy it carries. 

4. In this paper one of the motivations for studying the injectivity problem 
was the thermoacoustic tomography. One wonders then if considera- 
tions of 2D problems (rather than 3D ones) bear any relevance for 
TAT. In fact, they do. If either the scanned sample is very thin, or the 
transducers are coUimated in such a way that they register the signals 
only coming parallel to a given plane, one arrives to a 2D problem. 

5. Most of our results can be generalized to some Riemannian manifolds, 
in particular to the hyperbolic plane (where the analog of Theorem [7\ 
has not been proven yet). We plan to address these issues elsewhere. 
E. T. Quinto has recently announced a version of Theorem ^1 in the 
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case of distributions for the spherical transforms on real-analytic Rie- 
mannian manifolds with infinite injectivity radius and an analytic set 
S of centers 

6. A closer inspection of the results of Section |21 shows that most of them 
have their local versions, where it is not required that the whole trans- 
form Rs of a function vanishes, but rather only for radii up to a certain 
value. One can see an example of a local uniqueness theorem for the 
circular transform in P5|. We hope to address this issue elsewhere. 

7. As J. Boman notified us during the April 2004 AMS meeting in Lawrenceville, 
he jointly with J. Sjostrand, being unaware of our work, had recently 
independently obtained some results analogous to some of those pre- 
sented here (e.g., to Theorem ITHj) . 

8. We have not touched the problem of finding explicit inversion formulas 
for the circular transforms. Such formulas are known for the spherical, 
planar, and cylindrical sets of centers |3[in[IIllI3ll2HlEIlEZlEHlE3- 
They come in two kinds: the ones involving expansions into special 
functions, and the ones of backprojection type. Exact backprojection 
type formulas are known for the planar geometry ^01 1^ and recently 
for the spherical geometry in odd dimensions ^2] if the function to be 
reconstructed is supported inside the sphere of transducers. 

Another problem deserving attention is finding the ranges of transforms 
Rs- Such knowledge could be used, for instance, to replenish missing 
data. Some necessary range conditions have been recently obtained in 
[221 for spherical location of transducers. 

An important problem of reconstruction with incomplete data was 
treated in ^25,i40j based on an earlier work by E. T. Quinto in pij . 

9. An important integral geometric technique of the so called K-operator 
has been developed in I. Gelfand's school (e.g., [ISlITn]). It has been 
applied recently to the problems of the circular Radon transform (see 
[T7j , the last chapter of ^B] , and references therein) , albeit applicability 
of this method to the problems of the kind we consider in this paper is 
not completely clear yet. 
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